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Abstract 

Coherent pulse control for quantum memory is viable in the optical domain but nascent in mi- 
crowave quantum circuits. We show how to realize coherent storage and on-demand pulse retrieval 
entirely within a superconducting circuit by exploiting and extending existing electromagnetically 
induced transparency technology in superconducting quantum circuits. Our scheme employs a 
linear array of superconducting artificial atoms coupled to a microwave transmission line. 
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Circuit quantum electrodynamics has successfully realized numerous quantum optical 
phenomena in the microwave (/xw) domain [1] and is a promising architecture for quantum 
information technology |2j. Various superconducting, artificial atoms (SAAs) have been 
realized including Cooper pair boxes {sj, transmons Q|, flux qubits [5:] and fluxonium jfj]. 
These SAAs can be coupled to a /xw transmission line that serves as a 'quantum bus' for 
/xw photons Q. For efficient large-scale quantum information processing, synchronization 
of information flow between highly separated SAAs requires the /xw photons to be stored in 
quantum memories [8] coherently and retrieved on-demand. 

One strategy is to store photons in high quality superconducting microwave resonators [9|, 



but they yield low fidelities 



10|. Another strategy is to couple the superconducting circuit 



to external spin-i systems so that quantum information stored in the /xw field can be 



transferred for on-demand storage and retrieval [12] . Exciting progress has been made in 
this direction by coupling the superconducting circuit to nitrogen vacancies in diamond [kJ 
but the storage time is still limited to ~ 30 ns [jjj], and hybridizing with another quantum 
information medium adds extra complexity to the fabrication process. Pulse storage may 
be achievable by controlling the transition frequencies of SAAs [l5| but this scheme lacks 
coherent control for quantum memory. We propose to create a linear array of SAAs coupled 
to the transmission line, with this array serving as a pulse storage and on-demand retrieval 
medium. Our scheme is similar to quantum memory approaches in atomic systems [l6| and 
leverages off recent experimental demonstrations of superconducting-circuit electromagnet- 
ically induced transparency (EIT) [17]. The point of our scheme is to create /xw quantum 
memory in a superconducting circuit without the need to hybridize media, akin to the optical 
case with atomic quantum memory jg]. 

Although various choices of SAAs are available, we explore using fluxoniums |6] because 



they can serve as highly-coherent three-level SAAs [18] . Three-level fluxoniums are particu- 
larly amenable to /xw-field control via Autler-Townes splitting or EIT 19] . As illustrated in 
Fig. H^a), in our model, n SAAs are periodically spaced with inter-SAA distance / optimized 
to yield maximum storage efficiency The total medium length is d = (n — 1)1. Aperiodic 
spacings affect the position of the scattering resonances in the EIT spectra but do not alter 
the efficacy of our scheme. We shall see that interference due to inter-SAA scattering mod- 
ifies the transmission spectrum in a geometry-dependent way, but these complex spectral 
features are not an obstacle to realizing /xw pulse storage and retrieval. 
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FIG. 1: (a) Linear array of fluxonium (rectangles with dashed side) periodically spaced with 
separation / and coupled to (gray) transmission line, (b) Three energy levels of the fluxonium with 
on-resonant control driving of the |e) f-> |s) transition with frequency oo c and off-resonant driving 
of the |e) -B- |g) transition with frequency oj s and detuning 5. (c) Incoming and outgoing fields 
(indicated by arrow direction) for one SAA labeled by field amplitudes En and E' L for incoming 
fields and £"l and E'^ for outgoing, with L and R denoting 'left' and 'right', respectively. 

The energy levels |g), |s) and |e) and driving fields for each three-level SAA are depicted 
in Fig. UXb). Two driving fields are required to realize EIT. The control-field frequency u c 
is resonant with u cs (with u xy the frequency for transition |x) -H- |y) and x,y G {g, s, e}) 
and a signal field of frequency u s detuned by 5 = u s — u cg from u eg . When the control field 
is off, the SAA is opaque to the signal field; turning on the control field makes the SAA 



20] 



transparent. In addition, this three-level SAA can yield 'EIT with amplification' 

Previous studies focusing on EIT in superconducting circuits concentrate on the case 
of a single SAA and predict [2l| and demonstrate [17| changes to the transmission and 
reflection characteristics of the single SAA but do not consider phenomena such as slow 
light and properties such as optical depth that are crucial for building quantum memory. 
Here we consider a linear array of n SAAs to enable on-demand storage and retrieval of 
/j,w pulses. Before the fiw signal pulse arrives, a large-amplitude long-duration /iw control 
field coherently pumps the medium from |s) to |e) via the transmission line, and the field 
addresses the SAAs equally when they are relaxed to |g). The pumping results in coherence 
between the SAA energy levels, which suppresses the absorption of the small-amplitude 
signal pulse as well as slowing the signal propagation for pulse storage. 

When the signal pulse has entered the SAA medium, the control field is turned off to 
absorb the photons, which are converted to dark-state polaritons for storage. The control 
field is subsequently turned on to retrieve the signal pulse on-demand. This method of fiw 



pulse storage in a quantum circuit is inspired by the technique for storing an optical pulse 



in an EIT-enabled medium 



22|, which means the artificial atomic medium requires large 



optical depth for high storage efficiency [8|, and we show that our array of SAAs is capable 
of delivering the requisite optical depth. 

As the medium is one-dimensional (ID), the incoming signal fields encountering a SAA 
are scattered into transmitted and reflected components as depicted in Fig. die). For a 
strong non-depleting control field with constant Rabi frequency modulus \Q\ throughout the 
medium, the reflection coefficient r(S) for the signal field scattering from a SAA is 23] 



r + r r 

1 eg T i es -p - 1 - sg 



; e 

AS) 



r e = ^^,r s = ^, (i) 



with r xy the decay rate for the |x) -H- |y) transition. By treating the SAAs as point scatterers 
with linear susceptibilities, each SAA acts as a boundary for the transmission and reflection 
components of the signal field. Therefore, we apply the reflection-and-transmission transfer 
matrix method |24j to calculate the transmission through a sequence of SAAs. The transfer 
matrix method has been employed for studying analogs of EIT in plasmonic systems but 
without considering pulse storage 25j. For a single SAA, the induced transformation of the 
signal field is 

\E L (S) ) \r(S) 1 - 2r(S) j \e' l (S) ) 
We include phase shift ip = Iuj s /c, taking into account the free propagation of the signal field 
between the SAAs through the transmission line at speed c. The overall transformation by 
the linear array, denoted as matrix M, is given by 

' EM ) = M W ( E 'M (3, 

for the subscript 1 denoting the first SAA and n the last SAA of the ra-SAA array. 

As E in := Ei R is the only nonzero input field, E' nL = 0, only the upper- left element of 
the response matrix 

_ {A + + B)(A. + B) n - (A + - B){A. - BY 
11 2 n+1 e i ( n - 1 '><fB(r - l) n 

A± = ±e 2iv {l-2r) -1, 

B =v /( e 2^ - l)(e 2 ^(l - 2r) 2 - 1) (4) 
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is needed to obtain the output signal field E OVLt = E- m /Mn, with transmission T = |l/Mn| 2 . 

Let us now consider an example of how this system works by judiciously choosing fa- 
vorable yet realistic parameters. The transition frequency values are u eg /27i = 10.4 GHz 
and u es /27r = 6.99 GHz at flux bias $ = 0.32$ f° r "^o = h/2e the flux quantum; 
this flux bias yields long coherence times 18) . The |e) -H- |g) transition wavelength is 
A = 2nc/u eg = 11.6 mm and c = 1.2 x 10 8 m/s 26|. We treat decoherence as being due 



fundamentally to spontaneous emission |23[. Using Manucharyan et al.'s 



18) measured re- 



laxation parameter T sg 0.167 MHz, we calculate the fluxonium relaxation parameters to 



be r eg = 173r sg and T es = 40r sg 



29]. 



The field transmission in a translucent medium is related to the optical depth a by 
T = e~ a . For equally spaced SAAs, the length of the medium increases as n increases and 
Beer's Law holds if a oc n. For the choice of parameters and for two choices of inter-SAA 
separation Z, at A/4 and at one other value, we can see in Fig. [2]^a) that the optical depth 
at 5 = with control field off (Q = 0) is a linear function of n in accordance with Beer's 
Law; additional modulation due to inter-SAA scattering interference is negligible. For even 



a dozen fluxoniums t 
EIT of less than 25 



re optical depth is high compared to typical optical depths in optical- 



27) . From equation 03]), we find that the medium absorbs the signal 



field maximally when I — > (2m + l)A/4, where m is any non-negative integer. An exception 
to Beer's Law occurs when I — > m\/2, in which a — > 21n ^ 1+ ^~|q^°^ ■ 

Inter-SAA scattering interference occurs during transmission of a /iw field through a 
periodic array of SAAs, which may have complications that are absent in optical-EIT. We 
explore the effect that inter-SAA scattering interference has on EIT within the ID array 
of fluxoniums as depicted in Figs. |3]^a) and[3^b). Due to interference effects, in Fig. |3]^a) 
the transmission profiles with and without the control field on are markedly different from 
Fig. EJ^b). Fig. Mjo) is plotted using identical conditions as for Fig. |3]^a) except with a 
different inter-SAA spacing. In Fig. [3]^a) scattering creates asymmetry and oscillations but 
the EIT features remain intact. Fig. G^b) looks like textbook EIT, with 1 — T replacing 
absorption, which is realized by choosing some appropriate inter-SAA spacings. 

In our scheme, the /iw input signal field is a pulse, and its speed is quantified by the 
group velocity. The transmitted signal pulse acquires a phase factor e lkd from propagation 
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FIG. 2: Characteristics of a periodically spaced SAA array with n fluxoniums and inter-SAA 
separation I. (a) Optical depth avsn for / = A/4 = 2.90 mm (gray dots) fitted to a straight line 
a = 4.33n and for I = 0.74 mm (black dots) fitted to straight line a = 2.78n. For (b-c) the relevant 
parameters are / = 1.50 mm and the control-field Rabi frequency is = 218 MHz with numerical 
evaluation (dots) and the optical-EIT formulae (dashed), (b) Group velocity of a pulse through 
the array, (c) EIT window width w vs n. Analytical evaluation (solid) of w using Eq. ([6]) shows 
oscillatory behavior. 
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FIG. 3: Transmission through n = 100 fluxoniums with inter-SAA separation (a) I = 0.3 mm and 
(b) 1 = 1 mm for = (dashed) and with O = 309 MHz (solid). 



through the medium, and the dispersion relation and group velocity are 



k(u s ) = -arg 



/dfcV 1 
l r""\M 11 (u s )j ,Vg ~ \du s 



(5) 



s/ 5=0 



respectively. Based on T = Ts=oexp(—5 2 / w 2 ) [28( with Ts =0 w 1, which holds for small 5 
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(i.e. in the vicinity of the EIT window), the EIT window width w is given by 

d 2 lnT 



w 2 = -2 



d5 2 



(6) 

5=0/ 

Using Eqs. fl5]) and ([6]), we determine the inter-SAA scattering- included group velocity and 
EIT window width. Furthermore, we determine the group velocity and the width for the 
case of negligible inter-SAA scattering as in optical EIT. When EIT occurs, within the EIT 
window, the reflection of signal due to each SAA is small. Therefore, for the case of negligible 
inter-SAA scattering, the overall transmission coefficient is M^ 1 = (1 — r) n e^ n_1 - )</3 and the 
resultant transmission is T rs e -2nRe(r)_ g ence we have 



l^l 2 ,o, 2 /(r c + 2r s )|fi|2-4r3 



w w 4r eg/ 3v^' " N V 2 re g (4r e S r s + \n\*)* ■ (7) 

For x the linear susceptibility, the overall transmission coefficient is = 

e iui s Xd/(2c) e i(n-l)ip 28] SO 



1 2nr r 



- + 7- T^,J • (8) 



Equations ((Tj) and (jSJ) are the ID version of the EIT window width and group velocity 
expressions from optical-EIT theory [2^|, and we refer to these two equations as optical- 
EIT formulae. We evaluate the scattering- included v g and w and compare them against 
optical-EIT formulae in Figs.[2]^b-c), respectively, to see how closely the optical-EIT formulae 
apply in the case of a /iw pulse propagating through a fluxonium array. The agreement 
between the inter-SAA scattering-free theory and scattering-included theory shows that, 
with the chosen parameters, inter-SAA scattering interference has negligible effect on the 
group velocity. Furthermore, for the case of including scattering, w is larger with w oc -7= 
plus an additional modulation due to inter-SAA scattering, which appears as oscillations 
that dampen exponentially as n increases [30|. 

To show that storage of a /iw pulse in the fluxonium array is indeed possible using the 
EIT phenomenon, we consider an input pulse whose spectral profile is Gaussian with respect 
to detuning 5 with 5 = corresponding to input frequency u eg . The spectral function for 
the dimensionless input field is thus a Gaussian function E in (S) with standard deviation a 
and normalized field strength according to J \E- m (5)\ 2 d5 = 1. In our numerical evaluation, 
we determine the pulse storage efficiency 17(71, fl, I, a), which is defined as the percentage of 
the pulse trapped inside the medium as soon as the control field is turned off. 



We now employ the following conditions to evaluate the efficiency under reasonable ex- 
perimental conditions. The medium must be sufficiently transparent when the pulse enters, 
we set Q such that T V ^ Q = 0.99 at 5 = 0. We require that « g <c when the control field is 
on, which is guaranteed by choosing I such that v g ~ c/100 from Eq. ([H]). Not only does this 
condition ensure that the pulse is slowed inside the medium but also assures that the control 
field can be switched off to make the medium opaque before the pulse escapes. The pulse 
must fit inside the EIT spectral window so we choose a such that J \E out (S)\ 2 d5 ~ 0.98 as a 
benchmark. If the pulse is temporally longer than the available time window d/v g , only the 
central part of the pulse within this time window is stored and the pulse storage efficiency 
is reduced. 

Storage efficiency vs n is shown in Fig. H^a). The corresponding values of SAA sep- 
aration I, control-field Rabi frequency Q and the standard deviation of the pulse a from 
the constraining conditions are shown in Figs. H^b-d) respectively. Although some of the 
separation values in Fig. 11(b) do lead to distortion in the EIT spectra as in Fig. |3(a), the 
effect on the efficiency is negligible with the corresponding pulse widths. As the number of 
SAAs increases, the pulse storage efficiency increases as expected. For five fluxoniums, the 
efficiency is r] w 0.15, and for n = 100, the efficiency is rj w 0.72, which is comparable to 
the (retrieval) efficiency observed in optical-EIT pulse storage medium 27J. Asymptotically, 
r) approaches near unity as n becomes large. For example we evaluate the efficiency for 
n = 300 and find that rj w 0.91. Therefore, pulse storage is possible with a few fluxoniums 
and asymptotically efficient with increasing number of fluxoniums. At 0.32$o flux bias, the 
coherence time for fluxonium is about 1 /xs 18j so our proposed memory could store a pulse 
for approximately 1 /xs with negligible loss when EIT is off. Loss increases with the number 
of fluxoniums when EIT is on but strengthening Q ameliorates this loss, see Fig. 4(c). The 
pulse is retrieved after storage by turning on the control field and is optimally retrieved if 
the retrieval pulse is the time-reversal of the input pulse [271 ]. 

In conclusion, we have shown that coherent storage and on-demand pulse retrieval are 
feasible using existing superconducting artificial atom technology ]6] provided that an array 
of sufficiently large number of the artificial atoms is constructed as suggested. Our scheme 
obviates the need to hybridize the technology, e.g. by coupling to nitrogen-vacancy dia- 



mond defects 



131 ] . Hybridizing is a promising approach but introduces unwanted complexity 



and the storage time is limited by inhomogeneous broadening. Our scheme is surprisingly 
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FIG. 4: (a) Optimal storage efficiency for n fluxoniums. Solid line is an interpolated fit. (b) 
Corresponding fluxonium separation I. (c) Corresponding control-field Rabi frequency, (d) Corre- 
sponding standard deviation of the Gaussian spectral profile of the input pulse. 

efficient even for few artificial atoms, for example delivering 15% efficiency with just five 
fluxoniums and with approximately 1 /is storage time. This efficiency is dramatically im- 
proved with more fluxoniums and ultimately reaches near unit efficiency asymptotically. 
Our scattering-free ID optical-EIT equations provide a good guideline for how to construct 
quantum memory in superconducting circuits. Numerical evaluations are then necessary to 
assure that the theory gives quantitative agreement with experiments. Our results are im- 
portant for information synchronization and quantum memory in superconducting circuits, 
which are some of the most promising quantum computing technologies. 
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with A. Blais and J. Joo. 
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I. SUPPLEMENTARY MATERIAL 



Inter-SAA scattering interference occurs in a one-dimensional array of superconducting 
artificial atoms, which leads to oscillations in the electromagnetically induced transparency 
window width w as a function of the number of fluxoniums n. Here we derive an expression 
for w(n) that shows the w oc ^= behavior with an additional oscillatory modulation term. At 
the vicinity of zero signal field detuning, i.e. S ~ 0, the transmission T(5) = T(0)exp(— ^) 
with T(0) ~ 1. Hence, the window width is 

2 



w 2 



"dP" 



5=0 



where the transmission T 



2 



For r the reflection coefficient for a single fluxonium 



Mil 

and the phase shift due to free propagation between a pair of adjacent fluxoniums, 

3=1,2 

A + + B 

2 2e-^5 ' 

F - A - + B 
1 2e**(r-l)' 

A_ - 

Fo = 

2e^(r-l)' 

A± =±e 2 ^(l-2r) - 1, 
B =v /(e 2 ^ - l)(e 2 ^(l - 2r) 2 - 1). 

Note that the square root in the expression for B may take a plus or minus sign. Changing 
the sign leads to a swap of expression for G\ and G 2 as well as for F\ and F 2 ; thus, the 
expression for Mn is the same for either case. In the following, we obtain an expression for 
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d 2 lnr 



-InT =ln GjFJ 1 + c.c, 
3=1,2 

dlnT 1 ^ / „ dlnF dG, 



^ / _ dlnF, dGA 

;— 1 on / 



<W J2j=l,2 G 3 F j 3 , 2 

d 2 lnT 1 /dlnFA 2 , /„ d 2 lnF,- dG 3 - dlnFA d 2 G^ 



n 2 + Gi—r^r 1 + 2^^± n + ^ F™ 



d ^ 2 Ej=i )2 ^ 2 V V d6 J V d ^ 2 d5 d5 J d5 2 



J > \j=l,2 

We define 



a, — Gj 



2 



dV} dG.dlnF, 
° J d<5 2 + ^ d5 d5 ' 

_d 2 G 3 

dG, 



d5 2 G 1 F? + G 2 F 2 n \ G 1 F{ l + G 2 F? 

' / (G 1 Fl l + G 2 F 2 ")(a 1 F 1 " + a 2 F 2 ") - + &' 2 F 2 ") 2 ) n 2 



J d5 
such that 

d 2 lnT ( ai n 2 + M + ci)Fp + (a 2 n 2 + M + c 2 )F™ / (b[n + c[)F? + (6' 2 n + c 2 )F 2 ™ x 2 



+ c.c. 



(GiF" + G 2 F 2 ") 2 

+ ((GiF" + G 2 F 2 ")(6 1 F 1 " + 5 2 F 2 ") - 2(&' 1 F 1 " + b'.F^F? + c 2 F 2 ")) n 
+ (GiF" + G 2 F 2 n )(ciF™ + c 2 F 2 ") - (c^F™ + c 2 F 2 n ) 2 ) + c.c. 

= 2Re ( (GWTW ( (fl2Gl + aiG2 " 2b%) n2 

+ ((61G1 - 2b' 1 c[)F 2n + (b 2 G 1 + 6iG 2 - 2b[c' 2 - 2b' 2 c[) + (b 2 G 2 - 2b' 2 c' 2 ) F 2n ) n 
+ (dd - c'^F 2 " + (c 2 G! + Cl G 2 - 2c' lC2 ) + (c 2 G 2 - c 2 2 )F| n )) 

as FiF 2 = 1, a\G\ = b' 2 and a 2 G 2 = b' 2 2 . For the chosen parameters (see paper), at 
5 = 0, we have F 1 = 0.689 + i0.725, F 2 = 0.689 - i0.725, G x = 3.33 x 10~ 9 + i3.50 x 10" 9 , 
G 2 = 0.689 + i0.725 and F 1 -G 2 = -6.94 x 10~ 5 - i7.30 x 10~ 5 ; thus, |Gi||Fi| n < |G 2 ||F 2 | n 
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for all n. Also, a 2 G 1 + a ± G 2 - 2b[b 2 = 1.46 x 10~ 27 - i2.88 x 1(T 26 , hd - 2b' 1 c' 1 = 
3.32 x 10~ 35 - il.49 x lO" 35 , b 2 G x + b x G 2 - 2b[c' 2 - 2b' 2 c[ = -2.87 x 10~ 23 + i5.66 x 10~ 22 , 
b 2 G 2 - 2b' 2 c' 2 = -1.42 x 10~ 18 + i6.38 x 10~ 19 , ad - c[ 2 = -3.40 x 10~ 28 + i6.71 x lO" 27 , 
c 2 G , i+ciG 2 -2c / 1 c 2 = 7.01xl0- 20 -il.39xl0- 18 and c 2 G 2 -c' 2 = -7.01 x l(T 20 +il.39x 10~ 18 . 
Therefore, the n 2 term and the _F 2n terms are insignificant. The simplified expression for 
the window width in Hz is 



y/7.10 x 10" 19 n+ 1.39 x 10" 18 (1 - (0.9999) 2n cosl.62n) ' 
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